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We introduce a generalized joint density of states (GJDOS), which incorporates the coherent 
factor into the JDOS, to study quasiparticle interference (QPI) in superconductors. The intimate 
relation between the Fourier-transformed local density of states and GJDOS is revealed: they corre- 
spond respectively to the real and imaginary parts of a generalized impurity-response function, and 
particularly share the same angular factors and singular boundaries, as seen from our approximate 
analytic results for d-wave superconductors. Remarkably, our numerical GJDOS analysis agrees well 
with the QPI patten of d-wave cuprates and s±-wave iron-based superconductors. Moreover, we 
illustrate that the present GJDOS scenario can uncover the sign features of the superconducting gap 
and thus can be used to explore the possible pairing symmetry of the K a; Fe2-j / Se2 superconductors. 

PACS numbers: 



In superconductors, scattering of Bogoliubov quasi- 
particles by impurities can form quasiparticle interfer- 
ence (QPI) patterns which can be examined by the 
Fourier transformed scanning tunneling spectroscopy 
(FT-STS)^. For high-T c cuprates the characteristic 
wave vectors as well as their dispersions measured by 
the FT-STS experiments are interpreted quite well by 
the octet model^, which is based on the joint density 
of states (JDOS) argument. However, JDOS alone is in- 
capable of explaining the variation of the QPI patterns 
when the magnetic field is applied. Indeed it was ob- 
served that QPI peaks corresponding to certain octet 
wave vectors are suppressed or enhanced by the magnetic 
field depending on whether the sign of the SC gaps is re- 
versed or preserved^. The reason lies in that the probabil- 
ity of the elastic scattering of Bogoliubov quasiparticles 
on the constant contour of energy (CCE) is governed not 
only by their JDOS but also by the matrix element of 
the impurity Hamiltonian. Considering that Bogoliubov 
quasiparticles in superconductors carry proper coherent 
factors reflecting their particle-hole dualism, this scatter- 
ing matrix element will bear a particular combination of 
the coherent factors relying on the type of the scattered. 
Taking into account an important role of coherent fac- 
tors, the sign-sensitive effect of the QPI pattern on the 
magnetic field is naturally explained in the sense that 
the Andreev scattering of quasiparticles is enhanced by 
the pair-breaking impurities (mainly vortices) induced by 
the field. Notably, this sign-sensitive effect has also been 
applied to verify the s± pairing symmetry for iron-based 
superconducting (SC) material Fe(Se,Te)£. 

The insufficiency of the JDOS argument motivates us 



to propose the generalized JDOS (GJDOS) in this paper, 
which intrinsically incorporates the coherence factors. 
We find that the GJDOS and the Fourier-transformed 
local density of states (FT-LDOS) can be derived from 
the same quantity. The close connection between GJ- 
DOS and the Fourier-transformed local density of states 
(FT-LDOS) are established, which confirmed by our an- 
alytic results of ci-wave superconductors under the Dirac 
cone approximation. Our numerical results of GJDOS 
capture the essential structure of the QPI pattern and 
explain its evolution with the magnetic field for both d- 
wave cuprates and s±-wave iron-based superconductors. 
Finally, to explore the pairing symmetry of K x Fe2- y Se2, 
we evaluate GJDOS numerically and find that the re- 
sults are distinct for the theoretically-proposed nodeless 
d-, nodal extended s- and s ++ -wave SC states. Thus the 
pairing symmetry of K a; Fe2-j,Se2 may be identified by 
checking the QPI pattern under the magnetic field via 
the FT-STS experiments. 

We first develop the basic formalism for both the FT- 
LDOS and GJDOS. For different types of weak impurity 
scatterers, the FT-LDOS can be divided into the corre- 
sponding channels under the Born approximation 7 , 

^(p, W ) = -ilm[VG(k + p,w+)a l G(k, W +] 11 (1) 
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iO . at are Pauli matrices in Nambu 



space and correspond to magnetic (i — 0), pair-breaking 
(i = 1, 2)- and nonmagnetic (i = 3) impurity potentials. 
G(k,z) = (z — ^kO"3 — Atci) -1 is the Nambu-Gor'kov 
Green's function. Here the SC gap Ak is a real function 
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for high-T c superconductors, and therefore 52 is zero ac- 
cording to Eq. ^ for the case of i = 2. Thus, we only 
discuss the i = 0,1,3 cases hereafter. 

In order to deliver the concept of GJDOS and see how 
the FT-LDOS is related to GJDOS, we introduce a gen- 
eralized impurity response function Li 



E 

s=±,k 



G lQ (k + S p, w +K Q %(k,w) (2) 



where A a (k, uj) = — ImG a i(k, ui + )/ir with Ax and A 2 
correspond to the imaginary part of the normal and 
anomalous Green's functions, respectively. From Li we 
can derive two quantities, 



Si(p,oj) = ReLj(p,w) 



Ji(p,w) = --ImZrj(p,£j) 

7T 



(3) 



where Si(p,ui) is identical to Eq. (JTJ, while Ji(p,w) rep- 
resents our GJDOS. It is clearly shown in Eq. that 
the intrinsic relation between FT-LDOS and GJDOS lies 
in that they are just the real and imagine part of the 
quantity Li(p,ui + ) in analogy to, e.g., the dielectric re- 
sponse function. In addition to the autocorrelation of 
the normal spectral function Ai—, GJDOS also involves 
the autocorrelation of the anomalous spectral function 
A2 as well as the correlation between Ai and A2 . There- 
fore GJDOS extends the usual JDOS used in previous 
works 2 *^— , and is more transparent than the compli- 
cated J(q, E, E') used in Ref.[f|. 

Now we show how the coherent factors are naturally 
and appropriately embodied in GJDOS. From A% (k, uj) — 
u^S(lo - _E k ) + v^8(u) + Ey) and A 2 (k, oj) = u^Vk[8(oj - 
E k ) - 5{uj + Ek)}, where u£ = (1 + v j = 

(1 - W^k)/2, u k w k - A k /2£ k , E k = V^k + A L and 
according to Eq. ([2]) and ([3]), we obtain for ui > 0^, 

Ji(p, u) = 2J2 Ci(k + p, k)5(u - E k+p )S(Lo - Ek) 

k 

(4) 

where Ci(k' , k) is the coherent factor as given in Table U 
for three types of impurity. Also shown are the sign sen- 

TABLE I: Coherent factors for three types of impurity. 



impurity type 


coherent factor 


enhanced 


suppressed 


nonmagnetic 
i = 3 


M k /M k (u k /M k - V k tV k ) 


+- 


++ 


magnetic 
i = 


M k /M k (u k /M k + V k ,V k ) 


++ 


+- 


pair-breaking 
i = 1 


u k ,u k (u k iv k + « k iv) 


++ 


+- 



sitivities of the coherent factors for different quasiparticle 
scattering processes from k to k', where the sign of the 
superconducting gap is preserved or reversed indicated 
by H — h ( ) or H — ( — h). For the pair-breaking impu- 
rity case, the coherent factor is proportional to A k + A k < 



for £ k ~ £ k ' , which leads to a striking effect that the scat- 
tering amplitude will be zero if A k and A k / having the 
same magnitude but opposite sign. This sign sensitivity 
of the coherent factor for different type of impurity has 
been observed in the vortex region of (i-wave cuprates 4 - 
and has also been applied to verify^ the s±-wave pairing 
symmetry of iron-based superconductor Fe(Se,Te). 

To further investigate the relation between FT-LDOS 
and GJDOS, we turn to derive some analytic results of 
the complex response function Li{p, z) for <i-wave super- 
conductors under the Dirac cone approximation. Around 
the node we make the linear expansions with £ k sa uf&i, 
A k « -t>A&2, and E k w y/ (vFki) 2 + (wa^) 2 - For the 
intra- node scattering, we obtain from Eq. fl2Jl, 
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i = 3 
i = 
i = 1 



where 8 devotes the angle between p and the [110] di- 
rection, and z = 2u/Ep. From the above equation and 
Eq. ©, we can see that FT-LDOS and GJDOS, although 
not proportional to each other, share the same angular 
factors and singular boundaries, which holds true even in 
the inter-node scattering case. These common features 
can explain why GJDOS is appropriate to explain QPI 
in superconductors. We find that our results of Li(p, ui) 
agree with those obtained from A(p,w) in Ref. [5||9| for 
the nonmagnetic and magnetic cases. 

We also find that our analytic result for intra-node 
scattering in d-wave superconductors by magnetic impu- 
rities is very similar to that for scattering in surface of 
3D topological insulators 14 and intra-valley scattering in 
graphene 1 ^ by nonmagnetic impurities. The physics lies 
in the common Dirac cone structure and the <tq scat- 
tering channel. For illustration, we express for d-wave 
superconductors in Nambu space i? k = £ k <73 + A k <7i ~ 
vpkia^ — Vj\k2<7i = — [(up/ci , 0, V^k-z) x cr] • ri2, which 
is equivalent locally to an anisotropic Rashba spin-orbit- 
coupling Hamiltonian for the surface states of topological 
insulators. 

The GJDOS for d-wave superconductors are calcu- 
lated numerically according to Eq. Q and the results 
are shown in Fig. Q] q 2 , q 3 , q 6 , q 7 (q 1; q 4 , q 5 ) are 
relatively brighter (darker) spots for the nonmagnetic 
case contrary to the magnetic case as shown in Fig. [ljb) 
and (c), which is consistent to Table fl] For the pair- 
breaking impurity, the diminishing of the GJDOS peaks 
at positions q2,q3,q6jq7 and enhancement at positions 
qi, q4, qs agree with the magnetic field dependence of the 
QPI pattern observed in high-T c cuprates 4 -. 

We now turn to address QPI patterns in iron-based 
superconductor^ 1 ^— . Let us first consider where QPI 
peaks are expected to appear. It should be noted that 
wave vectors that connect maximum density of state on 
CCEs as in (i-wave cuprates are not expected here for the 
lack of strong anisotropy of CCEs for iron-based super- 
conductors. However, nearly nesting of CCEs are typical 
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FIG. 1: (Color online) (a) CCE in d-wave superconductors 
under Dirac cone approximation. The GJDOS patterns calcu- 
lated according to Eq. Q for (b) nonmagnetic (c) magnetic 
and (d) pair-breaking [i = 1) impurity cases. The normal- 
state dispersion is £k = — 2i (cos k x +COS k y ) —At' cos k x cos k y — 
fi and the gap function Ak = 2Ao(cosfc !/ — cosk x ) with pa- 
rameters t = l,t' = -0.3, fx = -1,A = 0.2, w = 0.4. 



S(x) 



_2a_ 



with 7/0 = 0.01 the broadening parameter. 



for iron-based superconductor and QPI peaks would ap- 
pear at those nesting vectors q where enormous scatter- 
ing processes k — > k + q on CCEs are involved. Those 
QPI peaks around nesting vectors are easily confused 
with peaks resulting from static orders such as the spin or 
charge density wave. However, QPI peaks are sensitive to 
coherent factors for different types of impurity depending 
on the relative gap sign for each pockets and can be sup- 
pressed or enhanced by introducing more pair-breaking 
scatterers after intentionally applying the magnetic field. 
This has been observed in the FT-STS experiment for 
Fe(Se,Te) that strong suppression of QPI peaks at (±7r, 0) 
and (0, ±7r) due to pair-breaking scattering occurs at the 
nesting vectors that connect the hole and electron pock- 
ets, supporting the s± paring symmetry proposed by the 
spin-fluctuation mechanism. In FigJ5]we show three QPI 
patterns for three types of impurity. And one can see 
that with increasing magnetic field and accordingly the 
strengthening of the pair-breaking scattering, the QPI 
pattern will change from FigJ^b) to Fig|3Jd), as observed 
by the experiment^. 

We then investigate the QPI pattern of the recently 
discovered K x Fe2- y Se2 superconductor. This iron-based 
superconductor is interesting because the hole pocket at 
the center is merged beneath the Fermi level and only the 
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FIG. 2: (Color online) (a) Fermi surface in Fe(Se,Te). The 
gray/white color represents the positive/negative sign of the 
SC gap with s±-wave symmetry. The GJDOS patterns for 
(b) nonmagnetic (c) magnetic and (d) pair-breaking (i — 1) 
impurity cases. The normal-state dispersion 213 is £k = £o — 
ii[(l — A) cos (k x — 7r) + (1 + A) cos k y ] — fj,i for k near (ir, 0), 
£k = eo — ti[(l + A) cos k x + (1 — A) cos(k y — ir)] — fii for k near 
(0, 7r), and £k = £o — £2 (cos k x + cos k y ) — fi2 for k near (0, 0). 
The gap function is A cos k x cos k y . The parameters are set as: 
£0 = 0.035, h = -0.68, %i = 0.85, m = -1.23, 1x2 = 1.54, A = 
0.3, A = 0.01, u = 0.01. 



electron pockets are left over. It calls for a re-evaluation 
of the spin fluctuation pairing mechanism due to the 
absence of the hole pocket and also raises a question 
about the possible superconducting pairing symmetry in 
K x Fe2- y Se2- Nodeless d 2 .2_„2-wave2i _ — , nodal extended 
s-wave22 and s++-wav o 21 ' 1 25 pairing symmetries have 
been suggested. While the last one is just a remain- 
der of s± with the same gap sign on the electron pock- 
ets, pairing states of the other two are fascinating for 
their complex gap structures. For the nodeless d-wave 
pairing, the gap function is slightly anisotropic without 
sign change around each pocket, while changes its sign 
between pockets, as schematically illustrated in upper- 
left panel of Fig. As for the extended s-wave pairing, 
the gap function around each pocket is similar to the d- 
wave pairing of the cuprates and this gap structure is 
rotated by 7r/2 between pockets, as shown in the middle- 
left panel of Fig. [3] Distinct gap structures of different 
pairing symmetries indicates featured QPI pattern for 
various types of impurity. In return, the pairing symme- 
try for K x Fe2-ySe2 can be deduced by checking the QPI 
patterns from FT-STS experiments. 

The normal-state dispersions of K :c Fe2-j,Se2 are de- 
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FIG. 3: (Color online) (a): Fermi surfaces in Ka;Fe2-i/Se2. The gray/white color represents positive/negative sign of SC gap 
with nodeless d- (1st row), extended s- (2nd row) and s++-wave (3rd row) pairing symmeties. The GJDOS pattern for (b) 
nonmagnetic (c) magnetic and (d) pair-breaking (i = 1) impurity cases. Parameters are set as: t = 0.7, n = —1.266, A = 
0.2, A = 0.01, 7? = 0.005, w = 0.01. 



scribed by £k = — t[(l — A) cos (k x — it) + (1 + A) cosfcy] — 
ji for the electron pocket located around (ir, 0) and 
£k = — t((l + A)cosfc K + (1 — A) cos (k y — n)) — fi for 
the electron pocket located around (0,7r). A accounts 
for the ellipticity of the pockets. Gap functions are 
set as: A(cos/c2; — cosk y ) for the nodeless d x 2_ y 2-wave, 
A (cos 

kx "i" cosfcy) for the nodal extended 5-wave, and 
A cos k x cos k y for the s++-wave pairing states. In Fig. [31 
we show the GJDOS patterns for three types of impu- 
rity with nodeless d-wave, extended s-wave and s++-wave 
pairing symmetries. While the three patterns seem to be 
similar for the nonmagnetic case, they are distinct for 
the pair-breaking case, especially around nesting vectors 
Q = (±7r, ±7r). Those changes can be summarized as 
follows: with stronger pair-breaking scattering (by in- 
creasing the magnetic field experimentally), the GJDOS 
magnitude at Q is strongly enhanced for s++-wave while 
deeply suppressed for the nodeless d-wave and extended 
s-wave pairing states due to the sign difference of SC 
gaps between the two electron pockets. Furthermore, 
the nodeless d can be differentiated from the extended 
s by examining the enhancement or suppression of QPI 
signal around the center. Thus by introducing more pair- 
breaking scattering, which can be realized by enhancing 
the magnetic field, one can deduce the pairing symme- 
try of K a; Fe2-j / Se2 by checking the magnetic-field depen- 



dence of the QPI patterns. 

In summary, we have introduced a new quantity GJ- 
DOS and employed it to study QPI in superconductors. 
It has been shown that the the GJDOS can faithfully 
reflect both the appropriate coherent factors and the sin- 
gular behavior of FT-LDOS. For d-wave superconductors 
with the Dirac cone approximation, analytic results have 
been obtained for both FT-LDOS and GJDOS, finding 
that they share the same angular factors (arising from 
the coherent factors) as well as the singular boundaries. 
This sharing, which may be rather general as we fur- 
ther argued, can account for why GJDOS is appropri- 
ate to explain QPI in superconductors. We have also 
shown that GJDOS is a powerful tool in revealing the 
SC gap structure, and our numeral results agree very 
well with the STS experiments for d-wave cuprates and 
s ± -wave Fe(Se,Te). We have further used the GJDOS 
approach to explore the possible pairing symmetry of 
K a: Fe2-j,Se2, suggesting that FT-STS experiments may 
be rather promising for revealing its pairing symmetry. 
Finally, we wish to pinpoint that the GJDOS approach 
may also be applied to QPI in topological insulator 2 ^, 
where the nontrivial spin-orbit coupling may bring in an 
interesting effective 'coherent factor'. 

The authors acknowledge a helpful communication 
with T. Hanaguri. The work was supported by the 



5 



RGC of Hong Kong under Grant No.HKU7055/09P and Science Foundation of China No.10674179. 
HKUST3/CRF/09, the URC fund of HKU, the Natural 



1 J. E. Hoffman, K. McElroy, D.-H. Lee, K. M Lang, 
H. Eisaki, S. Uchida, and J. C. Davis, Science 297, 1148 
(2002) 

2 K. McElroy, R. W. Simmonds, J. E. Hoffman, D.-H. Lee, 
J. Orenstein, H. Eisaki, S. Uchida and J. C. Davis, Nature 
(London) 422, 592 (2003) 

3 Q.-H. Wang and D.-H. Lee, Phys. Rev. B 67, 020511 (2003) 

4 T. Hanaguri, S. Niitaka, K. Kuroki and H. Takagi, Science 
328, 474 (2010) 

5 T. Pereg-Barnea and M. Franz, Phys. Rev. B 68, 
180506(R) (2003); ibid. 78, 020509(R) (2008) 

6 T. Hanaguri, Y. Kohsaka, M. Ono, M. Maltseva, P. Cole- 
man, I. Yamada, M. Azuma, M. Takano, K. Ohishi, and 
H. Takagi, Science, 323, 923 (2009) 

7 L. Capriotti, D. J. Scalapino and R. D. Sedgewick, 
Phys. Rev. B 68, 014508 (2003) 

8 K. McElroy, G.-H. Gweon, S. Y. Zhou, J. Graf, S. Uchida, 
H. Eisaki, H. Takagi, T. Sasagawa, D.-H. Lee and A. Lan- 
zara, Phys. Rev. Lett. 96, 067005 (2006) 

9 T. Pereg-Barnea and M. Franz, Int. J. Mod. Phys. B 19, 
731 (2005) 

10 Q. Han, D.-B. Zhang, Z. D. Wang, arXiv:1103.0"345l 

11 Here (Jo + Js)/2 corresponds just to the usual JDOS, 
namely the autocorrelation function of the single-particle 
spectral function^. 

12 M. Maltseva and P. Coleman, Phys. Rev. B 80, 144514 
(2009) 

13 It is readily to derive the corresponding formula for uj < 0. 

14 H.-M. Guo and M. Franz, Phys. Rev. B 81, 041102(R) 



(2010) 

15 T. Pereg-Barnea and A. H. MacDonald, Phys. Rev. B 78, 
014201 (2008) 

16 Y.-Y. Zhang, C. Fang, X. Zhou, K. Seo, W.-F. Tsai, 
B. A. Bernevig, and J. Hu, Phys. Rev. B 80, 094528 (2009) 

17 J. Knolle, I. Eremin, A. Akbari, and R. Moessner, 
Phys. Rev. Lett. 104, 257001 (2010) 

18 S. Sykora, and P. Coleman, Phys. Rev. B 84, 054501 (2011) 

19 T. Hanke, S. Sykora, R. Schlegel, D. Baumann, L. Har- 
nagea, S. Wurmehl, M. Daghofer, B. Biichner, J. van den 
Brink, and C. Hess, Phys. Rev. Lett. 108, 127001 (2012) 

20 A. Akbari, J. Knolle, I. Eremin, R. Moessner, Phys. Rev. B 
82, 224506 (2010) 

21 K. Seo, B. A. Bernevig and J. Hu, Phys. Rev. Lett. 101, 
206404 (2008) 

22 F. Wang, F. Yang, M. Gao, Z.-Y. Lu, T. Xiang and D.- 
H. Lee, Europhys. Lett. 93, 57003 (2011) 

23 T. A. Maier, S. Graser, P. J. Hirschfeld, and 
D. J. Scalapino, Phys. Rev. B 83, 100515(R) (2011) 

24 Y. Zhang, L. X. Yang, M. Xu, Z. R. Ye, F. Chen, C. He, 
H. C. Xu, J. Jiang, B. P. Xie, J. J. Ying, X. F. Wang, 
X. H. Chen, J. Hu, M. Matsunami, S. Kimura, D. L. Feng, 
Nature Materials 10, 273-277 (2011) 

25 T. Saito, S. Onari, and H. Kontani, Phys. Rev. B 83, 
140512(R) (2011) 

26 H. Beidenkopf, P. Roushan, J. Seo, L. Gorman, I. Drozdov, 
Y. S. Hor, R. J. Cava and A. Yazdani, Nature Physics 7, 
939-943 (2011) 



